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Abstract 



We consider the one-dimensional Dirac equation for the harmonic oscillator and the associ- 
ated second order separated operators giving the resonances of the problem by complex dilation. 

1 -Q ' The same operators have unique extensions as closed PT-symmctric operators defining infinite 

C^ , positive energy levels converging to the Schrodinger ones as c — ^ cxd. Such energy levels and 

their eigenfunctions give directly a definite choice of metastable states of the problem. Precise 
numerical computations shows that these levels coincide with the positions of the resonances 
up to the order of the width. Similar results are found for the Klein-Gordon oscillators, and in 
this case there is an infinite number of dynamics and the eigenvalues and eigenvectors of the 

fyr^ • _PT-symmetric operators give metastable states for each dynamics. 

Cn '. PACS: 03.65.Pm, 03.65. Ge 
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O ! 1 Introduction. 



The harmonic oscillator is one of the fundamental dynamical systems both in classical and quantum 
theory [6j and the natural relativistic extensions have been investigated since a long time. In 



^ • the relativistic quantum mechanics literature, therefore, one can find many different models that 



are considered relativistic oscillators [IHTj. These share the common property of reducing to the 
usual quantum oscillator in the Schrodinger limit but, otherwise, they present large differences 
in the way the interaction is constructed. Obviously different models have different dynamical 
and spectral properties: for instance scalar potentials have been added to vector potentials in a 
phenomenological way with the purpose of stabilizing the system and giving confinement, [HUllj. 
or else the interaction has been obtained by adding to the free Dirac Hamiltonian a linear term 
in the position coordinate twisted by the Dirac /3 matrix |12H18| . The relevance of such models 
in physics is clear and goes beyond their use for an elementary description of the confinement 
mechanism, in view of the connection with QFT and of the application of the Dirac equations to 
new materials [19], where the relativistic effects are more easily measured. 

We here call spinor and scalar relativistic oscillator, respectively, the Dirac and the Klein- 
Gordon equations with a quadratic electrostatic potential because these models are very realistic 
and, in principle, can be experimentally implemented. Although they do not satisfy the necessary 
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conditions for the exactness of the Foldy-Wouthuysen transformations [20112T], it is obvious that 
the weh known harmonic oscillator is reproduced in the Schrodinger limit. Despite their apparent 
simplicity, however, the study of the relativistic electrostatic oscillators is not so easy as it looks at 
first glance, but the results are satisfying. 

Common wisdom tells us that strong instability is induced in quantum relativistic equations 
with unbounded potentials in the way described by the Klein paradox, [22] . Because of such insta- 
bility the strongly bounded relativistic systems should present continuous spectrum and complex 
resonances. This picture is well signified by a result that Plesset established in 1932 [23j and that 
can be summarized, with a rather paradoxical formulation, as follows: the Dirac equation with a 
minimal coupling that involves an electrostatic "bounding" potential does not admit bound states. 
Of course, since bound states are present in the Schrodinger limit, the natural question that arises 
is what happens in the transition from a non-relativistic to a relativistic regime and in which sense 
this transition can be considered continuous in 1/c. A first answer was given by Titchmarsh in |23] 
who studied in detail the perturbation treatment of the Dirac equation with piecewise linear vector 
potential. Due to the very cumbersome analytical calculations involved, the treatment was essen- 
tially restricted to the first perturbation order giving, however some more explicit informations on 
the properties of the continuous spectrum by investigating the Weyl function m(A) of the singu- 
lar boundary value problem [25] in the complex plane. Very recently numerical investigations of 
the Dirac equation with a linear and a quadratic electrostatic potential were presented in |261I27| . 
looking for a dissipative model: metastable states were found and the Schwinger pair production 
rate ^28j was calculated in terms of the spectral concentration and in terms of the imaginary part of 
the resonances of the Dirac equation in external linear and quadratic potentials. In [26j the density 
of the states [29] was determined at finite values of (1/c), finding a sum of Breit-Wigner lines whose 
width reproduced the pair production rate. In [27] the Dirac equation with a quadratic potential 
was revisited in the framework of the large order perturbation theory. The singular problem gave 
rise to asymptotic non oscillating series and the sum was calculated by the Distributional Borel 
(DB) method, |3L)H34j . coupled with a Fade approximation for the DB transform, finding complete 
agreement with [2U]. Moreover the spectral problem was clearly formulated stressing the reasons 
for the dissipative choice. It was also shown that this choice leads to the study of non self-adjoint 
operators defined by complex dilation: the imaginary part of the complex energy levels, interpreted 
as resonances, again reproduces the pair production. 

Although exact bound states do not exist, in this paper we show that particular metastable 
states can be identified. They are uniquely obtained from eigenvalues of the decoupled equations 
of Titchmarsh. At least for a small ratio of the interaction to the mass energy they are well 
approximated. The role of the corresponding levels is similar to the role of the real part of the 
resonances. In this case we have both and it is our purpose to compare them. Actually, for 
the metastable states we have more exact results. In fact the Klein-Gordon like Hamiltonians, 
obtained by Titchmarsh in [25] (see also [35] ) from the separation of the Dirac system of equations 
- and hereby referred as to Titchmarsh Hamiltonians - are uniquely defined closed operators 
with a discrete spectrum that gives the metastable levels of the Dirac equation. The Titchmarsh 
Hamiltonians exhibit what in current literature is known as the PT-symmetry, where P is the 
parity transformation and T the complex conjugation [Tl[36l[37]: our paper, therefore, adds also a 
meaningful contribution to the discussion on the physical meaning of the PT-symmetric operators 
that still keep arising great interest (see, for instance, |38H42| ). The metastable levels are studied 
in the next two sections. In section 2 we use the methods of the functional analysis supported 
by some recent results coming from the study of anharmonic quantum oscillators |361I43| and we 



give formal statements of our analytical results. We also show the possible role of the nietastable 
levels as stationary levels of the separated PT-symmetric dynamics. In section 3 we then present 
the complex extensions of the position variables we have used, that have been introduced for two 
reasons: a better variational approximation and the connection with known results, mainly those 
of references |36U43j . In the final section we present a discussion of the numerical results connecting 
the metastable levels to the resonances of the model. We have thus calculated some perfectly 
defined energy levels of the Titchmarsh Hamiltonians stable at the Schrodinger limit and we have 
established their relationship to the real part of the resonances determined in [27J by the the DB 
sum: we have thus proved that the difference is of the second order in the imaginary part of the 
resonances themselves and it can therefore be ascribed to the pair production rate, as it had to be 
expected on a physical ground. We have found that the best numerical way, reaching the very high 
precision necessary in order to compare the asymptotic behavior in a parameter Q ~ 0{l/(?) — t- 
with the DB sum of the power series expansion [27] is provided by a specialization of the Rayleigh- 
Ritz scheme [H] obtained through the matrix moments method |451 l46j. This method could also 
be used for higher values of Jl where, however, other approaches could be equally or even more 
efficient [261117] 



2 The Dirac and the Klein- Gordon one-dimensional oscillators. 

Let us consider the one-dimensional Dirac equation in an electrostatic potential V{x) |24p35j. 
Using the notations of ref. pi] we assume a two-component spinor wave function of the form 
X = ^[X2{x) , Xi{x)\ so that the explicit form of the Dirac equation reads 

-iW + mc^-V{x))Xi{x)-h — X2{x) = 
c ax 

h^Xi{x) + -{W -mc^ -V{x))X2{x) = (2.1) 

ax c 

We therefore see that the 'large' component of the spinor is X2{x). In the following we will assume 

1 



V{x) = ^ muj^x^. 



We rescale the spatial coordinate as 



,muj. 1/2 
X I— ;■ ( ^— ) X 
n 



so to make it dimensionless and we define the equally dimensionless parameters 

Prom ()2.ip we then get the system 

^ X2{x) -Vn{E-x^ + ^) Xi{x) = 

^X^{x) + V^{E-x^)X2{x) = (2.2) 



In order to reduce the system (12. 2p to separate second order equations, it is convenient to 
introduce the hnear combinations 

i;+ix) = ^{X,{x) + iX2{x)) i;.ix) = -^{Xiix)-iX2{x)), (2.3) 



yielding 



d i 

— V'+lx) -iVn{E- x2)^+(x) - —^ (V+(x) + ii^-ix)) = 



— V'-(x) + iVn {E - x^)^_(x) + —^ (V'+(x) + ii^-{x)) = . (2.4) 

dx 2\'Q 



Solving the first equation (|2.4p in ip^{x), 

V'-(x) = -2Vn-^i;+(x) + i(2n(E -x^) + l) V'+(x) 
dx 

and substituting in the second equation, we finally find for ip+{x) the second order equation 

H+{n,E)i^+{x) = Xi;+{x) 
where A = A(J7, E) = E + Q,E^ and where the Titchmarsh Hamiltonian 

H+(n, E) = -—^ - 2iVnx + (1 + 2En)x^ - Ox^ (2.5) 

is a formal operator in the Hilbert space of the states Ti = L^(]R). Here and in the following we 
will always take vO > 0. 

The assumption of a real energy E will prove to be self-consistent with the boundary value 
problem that turns out to be complete, in the sense that the formal operator H^{i},E) extends 
uniquely to a closed Hamiltonian H-^- = H-^-{il.,E) (with discrete spectrum). The unstable term 
—Qx^ is still present in H^, suggesting that the Klein paradox arguments could be brought to 
bear: the imaginary term —2i\/Ux, however, makes the problem complete and the PT-symmetric 
extension of H^ unique. Moreover the unstable quartic term can be easily dealt with. Indeed 
in |36j a direct demonstration was given of the equality of the eigenvalues of a pair of quantum 
systems: the double well and the unstable anharmonic oscillator defined by complex translation, 
both of them, obviously, making sense and being well defined. In |43J, then, a general proof was 
presented that the spectrum of H^ is discrete and positive: 

a{H+) = {Xn > 0}neN ■ 

The equivalent adjoint Hamiltonian is 

H_ = PH+P = TH+T = Hi , (2.6) 

where P is the parity transform and T is the time reversal. 

Let us also notice that the non self-adjointness of the Hamiltonians has a physical meaning: it 
allows to define uniquely and to distinguish the H± operators. Finally, the labeling of the levels is 
the same we have at the Schrodinger limit, in the hypothesis that only one level is associated to 
one eigenvalue (see equation (|2.8|) here below). We shall argue that this is the case and, moreover. 



we expect and we prove the existence of infinite positive levels, stable at fi = 0. Let us now express 
our results in a formal way. 

Theorem 2.1 Let fi > 0, i? G R, then it is uniquely defined the closed PT -symmetric extension 
with discrete spectrum, Hj^{Vt,E) of the operator H-^-{0,,E). The closed operator H^(Q,E) is then 
defined from i\2.6\\ . 

Proof. We uniquely define the operator H^ by the L^ conditions at ±oo. The asymptotic 
behavior of the fundamental solution for x — )• +oo is 

^+(x) ~ — exp(^iS{x) — ln(x)) = —^ exp(i5(x)) , 



where 



S{x) = I' y'ny^-{l + 2En)y^ dy ^ VU^ - ^^^J^^^ x, (2.7) 

the other having the behavior 

^_(x) ~ — exp(iS'(x)) as x — )• — oo . □ 

Theorem 2.2 Let E £ M, Q > 0, and fix the index state it, there are infinite positive simple 
eigenvalues Xn{^,E) n £ N of the Titchm,arsh operators for the harmonic oscillator. For each 
r^ > 0, there are infinite positive energy levels En{0,) of the spectrum of the Titchmarsh operators. 

Proof. Let us fix 17 > 0, n G N and recall that, according to [l3], for any E G M we have 
positive eigenvalues Xn{E) of H±{E). We thus signify by "level" any generalized eigenvalue En > 0, 
solution of the implicit equation QE"^ + E = Xn{E), namely 

E = ^{Vl + 'iXn{E)n-l) (2.8) 

The existence of a solution of (j2.8p is proved in the following way. By real rescaling, and using the 
perturbation theory, we obtain, for fixed n, Q > and large positive E, 

\n{E) = Vl + 2E9. (2n + 1 + 0{Vl/{1 + 2E^f/'^ ) <^ E + QE^. (2.9) 

The same behavior applies for E > and small fi, so that in this regime the solution of (j2.8p is 
unique. □ 

Let us now consider how these results can be brought to bear to the study of the Dirac equa- 
tion. From ()2.3p we see that the components Xi{x) and X2{x) are linear combinations of two 
solutions ip±{x) of the second order equations. The eigenstates of the two Titchmarsh operators 
are independent because they are complex conjugate with non-vanishing real and imaginary parts, 
as it can be seen from the equations and their asymptotic behaviors. If we take a pair of eigen- 
functions {ip+{x),ip-{x)) of the separated equations such that ijj+ix) = iV'-(x), this pair gives 
directly metastable states of the full problem, whereas the resonances give metastable states only 
by a cut-off: it is relevant that we have more exact informations |36ll43j on the metastable states 
than on the resonances. In the non-relativistic limit — )• 0, the energy levels -E'n(^) tend to the 
eigenvalues En of the Schrodinger Hamiltonian and both the states tend to the real Schrodinger 
eigenstate ipn corresponding to En- Recalling again (j2.3p we also remark that, for the metastable 
states, X2 ~ Re V'+ and Xi ~ Im^+. The sizes of the large and small components of the spinors, 
therefore, turn out to be different by many orders of magnitude, providing an approximation to 



levels of the metastable states which is very natural and more accurate than the one we could 
obtain starting from the Schrodinger eigenvalues and perturbing them by means of the relativistic 
terms: indeed we the metastable levels are approximated by the DBS that gives the sum of the 
complete perturbation series in (1/c^) . 
Take then the real operator 

K{n,E)= H{n,E) + 2iVUx. (2.10) 

This is just the Klein-Gordon Hamiltonian with quadratic electrostatic potential and it is defined 
as a closed PT-symmetric operator K^{Q,E) by the behavior of the fundamental solutions, 

^±(2;) ~ TT expUSix)) (2.11) 

\x\ 

as X — >• ±00 . K^ has positive eigenvalues A„ for any real parameter E [36','43] . Again the equivalent 
adjoint Hamiltonian, 

K^ = PK+P = TK+T = R'l 

is obtained by parity transform. Following the same proof of Theorem 2, we can show that there 
are infinite positive energy levels as in the Dirac case. The boundary conditions are similar to the 
resonance ones, but there are always mixed conditions Gamow, anti Gamow at ±00. We define 
as the PT-symmetric pair of dynamics the dynamics generated by the two Hamiltonians K±. It 
is interesting to note that the energy levels are the same for the two dynamics, and that the two 
states are complex conjugate. 

Remark 2.1 Let us stress that the small difference between H^ and K^ is, however, relevant: 
the formal operator K is in fact not uniquely implemented as a closed operator. In particular, we 
have infinite self adjoint extensions with discrete spectrum. Moreover, there is the closed extension 
defined by the Gamow condition at ±00, with eigenvalues considered resonances, in a generalized 
sense, given by the Distributional Borel sum of the perturbation series. 

In perfect analogy with the H± Hamiltonians, we consider a pair of isospectral PT-symmetric 
Hamiltonians K± with positive simple eigenvalues. In this problem, where the physical dynamics 
does not exists because is not unique, it is appropriate to consider as physical levels, in first 
approximation, the values taken from the common eigenvalues of the pair of Hamiltonians K±. 
The corresponding eigenstates of the pair of Hamiltonians K± and their combinations, as the real 
means, are metastable states for all the infinite dynamics of the problem. The real part of the 
resonances looks less physical of the eigenvalues of K±, even if are very similar for small Q. Let us 
remark that probably it is not appropriate to impose the same principles of QM to the relativistic 
QM, where the pair production effect exists. But it is appropriate to pay the minimum price. 

For real energy E, we have computed the eigenvalues of ()2.5p by a variational method called 
the matrix moment method. Since each eigenvalue A^ is positive, [43j, we get the corresponding 
positive energy level ()2.8p as in the Dirac case. 



3 Definition of the Hamiltonians by subdominant behavior on two 
disjoint sectors. 

In this section we will give the prescription for a non-ambiguous definition of the possible operators 
we can obtain from the formal operator H{Q, E) with positive parameters (V", E). We follow the 



methods and the terminology of references [36^148) . Let us fix arg(ix) = (j), and consider the six 
sectors: 

Sj = {-7T/6 <{^- J7T/3) < 7r/6}, -2<j< 3. 

In particular, we define subdominant in the sector Sj, —2 < j < 3, the solution '^j of the second 
order equation, H{il.,E)^j = 0, with the principal behavior for cp = J7r/3, 

I |3 

ln{^j{x)) = -Vn^-^{l + 0{\x\-^))as |x| ^ cx). (3.1) 

In certain cases, the L^ behavior extends partially to the closure Sj of the angular sectors. In 
particular, the PT-symmetric Hamiltonian H^{Q,E) is defined by the subdominant behavior on 
the pair of sectors (S'_2,5'2)- Actually, the solutions respectively subdominant on the two sectors 
can be identified with the two fundamental solutions at ±00. 

Let us now examine the complex contours that can be taken according to the operators we want 
to define. 

(i) Complex translation and distortion. Let us start again from H = H^(Q,E) for positive 
{y/n,E), as above. Consider next the complex translation 

ijj{x) I— 7- Tzip{x) = '4}{x + z) , 

and the translated operator 

HHn, E) = — — r - Q.E'^ - E- 2iVTi(x + z) + (1 + 2En)(x + zf - n(x + z)^ . (3.2) 
For z = iy, y > 0, H^ is uniquely defined as a closed operator by the fundamental solutions 

^^(x)~^ exp(i5(x + z)), (3.3) 

as X —7- ±00, respectively. A relevant observation is in order. As the fundamental solutions are 
coincident with the subdominant solutions on the pair of sectors (5-2, 5*2), we can prove that all 
the translated operators are isospectral with the operator defined by the subdominant condition 
on the two sectors above. Indeed it is easy to see that the fundamental solutions of the translated 
operators are the translated of the fundamental solutions of if+(/3,-E), with the behavior 

^±{x)r^-^e^p{iS{x)), (3.4) 

as X — )• ±00, respectively. Since the eigenvalues are the zeros of the Wronskian of the two funda- 
mental solutions, we can conclude that the translated operators are isospectral with the operator 
H^{Q,E) itself. In the numerical calculations the translation parameter y will be used as a vari- 
ational parameter in the application of the Rayleigh-Ritz method. We could also consider the 
PT-symmetry conserving distortion of entire functions, 

'(/'(x) I— ;■ Dip^x) = \/z'{x) ■(/'(z(x)), 

where z{x) = x(l + i tan(7r/6) x Vl + x^ ) , [49], yielding strictly sectorial operators, 

H^{VL,E) = DH+{n,E)D~^. 



These will be treated elsewhere. 

(a) Complex dilation. Let H = H-^-{Q,E) as above, with ^2 > 0. Consider first the complex 
dilation 



■ip{x) I— ;■ Tg^{x) = ^(xexp(i^)) = ip (x), 



and the operator 



He{n,E) = TeH{n,E)T^e, (3.5) 

where ^(x) is restricted to the dense set of L^ functions analytic on an angular sector 

< tan(6'o) with < \e\ < Oq < vr/G . 



Rez 

Such operators are defined by the subdominant behavior on the two pairs of sectors (5_2, Si) and 
{S-i, S2) and they are isospectral to the dilated ones with negative and positive 9 respectively. For 
positive 9 the fundamental solutions at ibc« have the behavior 

ln($^(2;)) = iS{xexp{i9)) + 0(ln(|2;|)), (3.6) 

as X ^ ±00, respectively. For negative 9 the fundamental solutions at ±00 are 

ln(*i(x)) = -iS{xexp{i9)) + 0(ln(|a;|)). (3.7) 



The two DB complex conjugate sums |30H34] apply to the energy levels of these operators defined 
by the subdominant behaviors on the two pairs of sectors (5_2,5'i) and {8-1,82), respectively. 
The application of the perturbation theory leading to the DB sums has been thoroughly studied 
in [27] by numerical methods, specializing, in particular, to the calculation of the pair production 
rate obtained from the imaginary part of the resonance. 



4 Results and conclusions. 

In order to compare the spectra of different boundary value problems, instead of integrating the 
differential equations as in [26j, we have calculated the eigenvalues using a slight modification of 
the classical Weinstein-Aronszajn method ^^ obtained from the theory of matrix moments, |45j . 
joint to the Newton procedure for finding the zeros of a function. In fact, due to the Jacobi form 
(j4.ip of the Hamiltonian matrix in the base of the occupation number states, this is equivalent to 
a Rayleigh-Ritz variational scheme made more efficient through the matrix moments that shorten 
the computational time and increase the precision of the calculation. 

We report, in Figure 1 the lowest energy levels of the Dirac electrostatic oscillator for different 
values of the parameter 0. 

We consider the translated, real and dilated operators uniquely defined on L^(M). These 
are compressed on the space Tin spanned by the first n Hermite eigenfunctions {V'jjjXn of the 
Schrodinger Hamiltonian Hq = p^ + (T^x^, where o" is a variational parameter. The matrix elements 
Hiii. of the fourth degree polynomial in x are obviously dealt with by using the recurrence relation 

xhk = (1/2) hk+i + k hk-i 
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Figure 1: The lowest energy levels of the Dirac oscillators for different values of il. Circles refer to 
the lowest positive energy level. Down triangles, squares and up triangles refer to the first, second 
and third excited levels respectively. 



for the Hermite polynomials hr = hr{x). The explicit relevant relations we are using read 

x^ hk = il/i) hk+2 + ik+l/2) hk + n{k-l) hk-2 

x^ hk = il/8) /ifc+3 + (3/4)(A:+l) /ifc+i + (3/2) k^ h^-i + k{k-l){k-2) hks 
x^ hk = (1/16) /ifc+4 + (l/4)(3+2A;) hk+2 + {'i/^){2k^ +2k+l) hk + k{2k-l){k-l) hk-2 
+ k{k-l){k-2){k-3) hk-i 

Thus the general form of the matrix reads 



(H., 



ik 



( Ao 


Bo 








... \ 


Co 


Ai 


Bi 










Ci 


A2 


B2 






(4.1) 



where the symbols An, Bn and C„ represent 4x4 blocks given by 



\An)ik 
{Bn)ik 
{Cn)ik 



K 



K 



4n+i,4n+k 
4n+i,4(ra+l)+fc 



H. 



4{n+l)+i,4n+fc 



K 



4n+fe,4(n+l)+i 



{Bn, 



ki 



with n = 0, 1, 2, . . . and i,k = 0, 1, 2, 3. 

The first 4n eigenvalues can be approximated by the eigenvalues of the Hamiltonian matrix 
truncated at the 4n-th order. From the matrix moment theory [l5] these eigenvalues are given by 
the zeroes of the determinant of the matrix polynomials P„ = Pni^), recursively defined in the 
following way: 



Po = h 



A. 



n-1 



Pn^l) , 



(4.2) 



I4 being the 4x4 identity matrix. The solutions of detP„(A) = are indeed the 4n-th order 
Rayleigh-Ritz approximants that give upper bounds for the eigenvalues. Although lower bounds 
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Figure 2: The plot of A = — ln(£^to ~ ^^{Edo)) vs. the number n of iterations of (|4.2p . Different 
symbols correspond to different values of Q. The horizontal dashed lines give 2Im(£'rfo) for the 
corresponding fi. 



for the eigenvalues do not prove to be useful in our present context, it could however be observed 
that they can be obtained by this same method, looking for the zeroes of the determinant of the 
matrix polynomial [46] 



P„(A)-P„(0)P-_\(0)P„_1(A) 

Let us now present the numerical results. We denote by "r, t, d" the quantities respectively 
related to the spectrum of the real, translated and dilated operators (|2.10p . (j3.2p and (|3.5p previ- 
ously discussed. In Table 1, for different values of fi, we give the complex d-eigenvalue Edo with 
lowest positive real part, reducing to the ground level in the non-relativistic limit. Since the effect 





n 




Eda 











0020 




0005017620+ i 


1 


17374083059 e- 


-144 





0025 




0006277579 +i 
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42079110945e- 
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0007539782 +i 
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72376665081 e- 


-96 





0035 




0008804241 +i 
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77543924661e- 


-83 





0040 
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00211928567e- 


-72 





0045 




0011339978+ i 
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08165603853e- 


-64 





0050 




0012611278+i 


5 


36447802132e- 


-58 



Table 1: The lowest d-eigenvalue E^q for varying fi. 

to be highlighted is really tiny, these eigenvalues have been calculated with great accuracy, both 
in the arithmetic precision and in the number of iteration of the recurrence relation (j4.2p . The 
largest value of the latter has been taken to be 300: although a large number of decimal figures is 
already stabilized by few iterations, a very high precision is however necessary for a comparison of 
the t with the d eigenvalues, that we show here below. As expected, the results obtained in [261127] 
when studying the resonances of the Dirac equation by the spectral concentration and by the DB 
sum, as well as the asymptotic behavior of the imaginary part, are confirmed by the much more 
precise data given in Table 1. In the application of the Rayleigh-Ritz method we have used as 
variational parameters the size of the imaginary translation y introduced in item (jS.Sp and the the 
'frequency' a of the operator Hq. There has been numerical evidence that the optimal values of 



10 



these parameters are 1 ^ y ^ 5 and 1 < cr < 2. In Figure 2 we plot 

A = -ln{Eto-Re{Edo)) 



for a different number n of iterations of ()4.2p in order to test the stabihzation of the data with n. 
We see that the saturation for decreasing values of fi requires increasing values of n: it appears, 
however, that n = 300 is already sufficient for 0, = 0.0030 and higher. The final value is almost 
coincident with minus twice the logarithm of the imaginary part of the eigenvalue Ed, represented 
by the dashed horizontal lines. 

Remark 4.1 This means that each energy level is given by the distributional Borel sum of the 
perturbation series modulo a correction of the second order on the pair production effect (see [27]). 

This property has been numerically checked by calculating the ratio A/(— 21nIm(ii^rfo)) for 
0.0030 < < 0.0050. A simple quartic Lagrangian interpolation on the data thus obtained gives 
0.999885 as a limiting value of the ratio. 

We finally want to consider the influence of the imaginary term proportional to vQ in (j3.2p . 
We have therefore calculated the difference of the first t and r eigenvalues and the ratio 

K={Etoin)-Erom)/^- 

A least square computation on the data for 0.002 < il < 0.005 gives k = 0.9999539755 + 
0.0284823904 , providing a numerical evidence that the difference of the previous eigenvalues 
vanishes in the limit of a vanishing Q. 

We have compared the differences of the first excited minus the fundamental level for the t and 
r operators: 

6 = [{Enin) - Etoin)) - (^.i(0) - E,,o(0))]/f). 

finding a more than linear vanishing behavior with vanishing fi. 

To conclude, in this paper we have completed the rigorous analysis of the basic model of 
the harmonic oscillator in ID Dirac and Klein-Gordon theory, using recent mathematical results 
concerning the anharmonic oscillators and considering the relevance of the physical symmetries 
[361137]. As we said, the result is satisfying, since it brings it one step closer to the QFT symmetries. 
We have also shown that discrete energy values, the exceptional points, can be obtained by the 
approximation of decoupling two equations, but without the presence of a scalar potential, whose 
physical interpretation is usually uncertain. The philosophy could be that: in relativistic QM it is 
not possible or convenient to keep all the principles of QM, but we can try to minimize the effects 
of the pair production. We have PT-symmetric energy levels comparable to the non relativistic 
ones. This way, the Klein Gordon equation appears more similar to the Dirac one. In this case, the 
PT-symmetric eigenvalues appears more physical than the infinite ones of the infinite selfadjoint 
Hamiltonians. The levels of the 2D and 3D Dirac oscillators will appear in a following paper. 
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